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We show that technique of Dyson equation in wave multiple scattering by spatially disordered 
discrete medium statistical theory leads directly to a dielectric permittivity tensor, which is char- 
acterized by spatial dispersion and obeys the generalized Lorentz-Lorenz formula. Introduced 
via this spatial dispersion an effective magnetic permeability demonstrates the diamagnetic 
property in limit of independent strongly reflected non-magnetic small spherical particles in 
accordance with earlier intuitive predictions. The revealed physical nature of the effective dia- 
magnetic property consists in that electric and magnetic dipoles induced in a particle by wave 
1-^ \ scattering give different contribution into the transverse and longitudinal components of the 

effective dielectric permittivity. Besides as appeared the diamagnetism under study is enhanced 
by appearance of additional effective dielectric displacement current in the medium. 

'■ 1 Introduction 

^ ; Effective medium theory or homogenization approach is a powerful way to handle the radiative 
properties of composite materials, such as, e.g., cosmic dust, aerosols, porous media, and mi- 
crostructured materials or "metamaterials". Several interesting potentials of these metamaterials 
have been demonstrated and/or suggested, including the realization of media with negative index 
of refraction [H[2] and artificial magnetism from conductors [3JHJE]. Even though some fundamen- 
tal concepts about the metamaterials such as negative permittivity and negative permeability are 
widely accepted and relatively well understood today, the problem to define unambiguously these 
parameters by systematic and clear way is not yet totally resolved. For example, the introducing 
artificial magnetism by [3j HJ E] is rather intuitive one. The two variable asymptotic expansion pro- 
cedure [6], which makes homogenization of Maxwell's equations in a periodically composite under 
condition that a period tends to zero, is presented as purely mathematic one and does not discuss 
the problem of artificial magnetism from non-magnetic conductors. The author [7] has derived 
a productive transformed integral equation for the electric field in periodically microstructured 
composite, with inhomogeneous term being proportional to electric field averaged over unit cell of 
the microstructure and integrating over unit cell in the integral term. This equation was solved 
numerically to get the spatially dispersive effective dielectric permittivity tensor and to consider, in 
particular, effective magnetic permeability via the spatial dispersion of dielectric permittivity ac- 
cording to [Sj- In [9], a semi-intuitive approach is constructed to solve analytically the transformed 
integral equation [7J, with obtaining some generalization of classical Lorentz-Lorenz formula for 
local electric field and similar formula for local magnetic field that gave possibility to obtain a 
corresponding generalization of the Clausius-Mossotti formulas for the classical "local" effective 
permittivity and permeability. 
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So far the approaches were discussed to homogenizing the Maxwell's equations in a periodically 
composites, with averaging over crystal cell. The aim of this paper is to show the systematic and 
clear way for effective magnetism with diamagnetic property definition in framework of electromag- 
netic wave multiple scattering theory by spatially disordered composites consisted of non-magnetic 
particles, each of which is characterized by given dielectric permittivity and conductivity. The wave 
multiple scattering theory [TUl [TTj starts with considering the wave scattering by a single particle 
that can be described in terms of particle tensor electric field scattering operator [12]. The prin- 
cipal point for our consideration is that this scattering operator is restored from Mie solution [13] 
for electromagnetic wave scattering field in all points of spherical particle outside, including near 
zone, in part of electric dipole and magnetic dipole scattering, that is with accounting the spheri- 
cal harmonics of unit angular momentum. As this takes place, the contribution of electric dipole 
scattering into electric field scattering operator in momentum space representation is appeared to 
be not dependent on momentums whereas the contribution of magnetic dipole scattering is a bi- 
linear function of scattering and incident momentums what is very suitable for effective magnetic 
permeability evaluation via effective dielectric permittivity tensor with quadratic spatial disper- 
sion [8] . The effective dielectric permittivity tensor is evaluated in our formalism from generalized 
Lorentz-Lorenz formula, which is derived by homogenization of stochastic Maxwell's equations for 
electromagnetic field in random medium under consideration with using the technique of Dyson 
equation pJJ[ [11] and averaging procedure over statistical ensemble of particles. Our observation 
the different contributions of a particle electric and magnetic dipole scattering into the transverse 
and longitudinal components of the effective dielectric permittivity is fundamental for appearance 
of an artificial effective diamagnetic properties in statistical ensemble of non-magnetic conducting 
particles. The authors [14] discuss similar phenomena in contribution of the electric-dipole and the 
magnetic-dipole transitions to the dielectric tensor in the framework of molecular light scattering. 

In what follows we present the materials of our paper making accent on technique of generalized 
Lorentz-Lorenz formula related to statistical ensemble of non-magnetic conductive particles that 
can be especially useful at appearance a method [15] of artificial dielectric disordered material 
manufacturing with controlled parameters. 



2 Homogenization of stochastic Maxwell's equations 

We start with stochastic Maxwell's equations for monochromatic electromagnetic wave field 

rotE = i-B, rotB = -i-£E (1) 
c c 

in homogeneous background (host) medium with frequency u, constant dielectric permittivity Eq 
and unit magnetic permeability /i = 1, where there are placed randomly particles with given dielec- 
tric permittivity £i, specific conductivity o\ and unit magnetic permeability. Symbol i(r) denotes 
stochastic complex dielectric permittivity of the medium under study, being equal to complex di- 
electric permittivity e± — E\ + Uttui/uj inside a single particle and e outside the particles. The 
problem of the random medium effective parameters is resolved by averaging the Eqs. ([TJ over par- 
ticle spatial distribution statistical ensemble, with the averaging denoted as (. . . ). In particular, 
an effective dielectric permittivity of the medium e e ff is obtained according to [TO], [Til HH1 EG] via 
a "decoupling" 

(*(r)E B (r)) = f dv>e%{v,v>) (E,(r')> (2) 

the averaged product (iE a ) of the random medium complex dielectric permittivity and the random 
electric field vector E a (r) with components a = 1, 2, 3. As one can see, effective dielectric per- 
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mittivity of the random medium is a tensor operator, with tensor indices being denoted by Greek 
subscripts and implied agreement about summing over repeated subscripts. The decoupling (j2J) is 
realized with the aid of the Dyson equation [TU1 E] , which is obtained via ensemble averaging the 
integral equation for stochastic electric field and written in a symbolic operator form as 

(E) = E {0) + G {0) M (E) (3) 

In the right hand side of this equation Ea (r) denotes the incident electromagnetic wave electric 
field; the tensor Green function of electric field in the host medium is given by G^(r) = (5 a p + 

y a V/3 Ao) G (r), with k = e l J 2 uj/c being the wave number in the host medium and G (r) = 
exp (ik r)/(— Anr) a scalar Green function; M a p(r,r') is the tensor mass operator. The effective 
dielectric permittivity and mass operator are connected between them by operator relation, e e ff/e 
= I - M/ky, where I is the identical operator. 



3 Lorentz— Lorenz formula 

The tensor Green function of electric field in the host homogeneous medium is strongly singular 
in the origin point and decomposed [TBI ITS] IT9] into delta Dirac function term 1/(3^)1 and the 
principal part G^°\ Substituting this singular tensor Green function into Dyson Eq.Q leads to a 

r 1 _1 

transformed mass operator that is defined by M = I — M and gives possibility to get for 

the effective dielectric permittivity an operator relation 

{e eff + 2e ) {e eff - e ) = -— = — u eff (4) 

which we call a generalized Lorentz-Lorenz formula. The transformed mass operator M is evalu- 
ated by applying the Feynman diagrams technique |T0l [11] or via equivalent method of asymptotic 
expansions [161 E] • The mass operator asymptotic expansions are obtained from asymptotic ex- 
pansions for ensemble averaged (T) electric field random scattering operator T of the medium 
and the equation (T) = M + MG^ (T). According to definition [12], the scattering operator 
T a p(r,r') satisfies the Lippmann-Schwinger type equation, T = V + VG^T, where the medium 
transformed scattering potential V is equal to sum of single particle transformed potentials Vi(r) 
= -&q 3 [ii(r) — sq] I [£i(r) + 2e ], provided that the particles are not overlapped. 



4 Generalized and classical Lorentz— Lorenz formula 

Making comparison between the obtained Lorentz-Lorenz formula (j3J) and the classical one [T3] 
shows that a quantity u e ff denoting the (j3J) right hand side has sense of the medium effective electric 
susceptibility per unit volume tensor operator. The ratio v e ff/fi, with f% being the particle number 
density, is an electric susceptibility per a scatterer. The left hand side of the generalized Lorentz- 
Lorenz formula is really looking like very familiar algebraic construction [13] and in the right 
hand side of this one stands the transformed mass operator, which includes all significant items 
influencing the dielectric and conducting properties of composites. Besides the frequency effects, i.e. 
accounting for intrinsic properties of constituency by their complex dielectric functions, the mass 
operator includes geometric effects, caused by resonant behavior of composite due to the shape and 
size of particles, as well as cluster effects because of correlations between particles, which can result 
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in significant modification of particle scattering properties. The obtained Lorentz-Lorenz formula 
includes such known cases as Maxwell-Garnett mixing rule with Clausius-Mossotti formula in the 
presence of multiple scattering on small correlated [18] and non-correlated [19] dielectric spheres. 
The case of non-correlated dielectric particles is described by transformed mass operator in the 
asymptotic form M ~ j dlfiTi where f\ is related to a particle center r x and T\ is the single 
particle electric field scattering operator with transformed scattering potential Vi(r — ri). 



5 Effective magnetic permeability 

If statistical ensemble of particles is homogeneous and isotropic in average, the operators e e ff 
and M are functions of spatial coordinates differences and can be expanded into Fourier integrals 
along spatial harmonics exp(ikr). Corresponding tensor Fourier transformations are defined by 
transversal t and longitudinal £ with respect to the wave vector k components, as E e ff )a p(\c) = 
(5 a/ 3 — k a k/3) e\ff(lz) + k a k/3 e e e jf(k), with k being the unit vector along k. In terms of these 
components the Lorentz-Lorenz formula (03) reads 



ef ff (k)-e M^{k) 



e¥ ff (k) + 2e 3fc 2 



(5) 



Actually, one does see in the left hand side of this formula the transversal and longitudinal effective 
dielectric permittivities [HI EH] for a monochromatic plane wave. 

The spatial transversal and longitudinal Fourier components of the transformed mass operator 
in the right hand side of Eq.(jS]) take in the case of non-correlated (independent) particles the 
form M t,e (k) = f\ T'^(k, k). Note, a quantity T Q( g(k, k') denotes the electric field tensor scattering 
operator of a single particle in spatial Fourier transform (momentum) representation with arbitrary 
wave vectors k and k'. After having been restricted on wave vector spherical surface k = k' = 
ko , the double transversal part of this tensor with respect to directions of both scattering and 
incident wave vectors gives the tensor scattering amplitude in far wave zone of the particle [12J. 
Nevertheless, we have met the tensor scattering operator T a p(k, k) for forward scattering direction, 
when the both wave vectors are identical and not equal to the wave number in the host medium with 
respect to their absolute magnitudes. We made supposition that this value of the tensor scattering 
operator is defined by its transversal and longitudinal components T*' (k, k) that is correct, for 
example, in the case of a spherical particle. Following [8] we denote 



x 1 _^ ei ff (k)-ej ff (k) 

Heff C 2 k 2 





M* (k) - M l (k) 


/k 2 


1 + M* (k) /3k 2 




1 + M l (k) /3k 2 



(6) 



setting next 47rl(k) = (1 — l//i e jj)(B(k)) and rewriting the ensemble averaged second Maxwell's 
Eq.([T]) identically as 

i k x «B(k)> - 47rl(k)) = --4ff(k) (E(k)> (7) 

c 

Eq.([6]) defines in the limit k — y an effective magnetic permeability [i e ff, with I(k) being a vector of 
effective medium magnetization. In result, Eq.^7]) takes a form of the macroscopic Maxwell's equa- 
tion with magnetic flux vector [20J, given by (H(k)) = (B(k)) - 47iT(k) and dielectric displacement 
current defined by longitudinal component of effective dielectric permittivity. 
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6 Electric dipole and magnetic dipole scattering 



Let us evaluate the effective magnetic permeability defined by (Q in the case of independent spher- 
ical particles. Instead to find a solution to the Lippmann-Schwinger equation we will construct 
here the single particle scattering operator T\ in approximation of electric dipole and magnetic 
dipole scattering. The incident electromagnetic wave induces into the particle the electric dipole 
and magnetic dipole moments p = eo^Eo and m = xH with dielectric and magnetic suscepti- 
bilities r] and Xi respectively. The Hertz's vector of electric dipole and magnetic dipole scattering 
Ili e ^(r) and Ua 9 \r) is given according to [21] by n« (r) = - (47r / so)p a Go(r) and ni ma ^(r) = - 

1 /2 

(47i/ikoe )e , i s 7 m / 3V 7 G'o( r ), respectively, where e a/ 3 7 is absolute anti-symmetric unit tensor. Bear- 
ing in mind the rule to evaluate the electric wave field via the Hertz' vector and also via the 
scattering operator, E = G^TiE^°\ one can get the following expressions for scattering operators 
of electric dipole and magnetic dipole scattering, respectively, in the form 

Ti7 9) (r,r') = -4n X e a ^e Kl V v V J 6(r)5(r'). (8) 



The evident expressions of spherical particle electric and magnetic susceptibilities rj and x are 
written out by comparison with Mie solution JT3J for electromagnetic wave scattering in all points 
of particle outside and have in Hulst [22] partial scattering amplitudes the form rj = (3i/2/cQ)eti and 
X = (3i/2kQ)bi. The transversal and longitudinal components of the obtained scattering operators 
in momentum representation for forward scattering direction are given by 

T (d) '(k,k) = T (eI)£ (k,k) = -4irr]k 2 ; 

T (mag)t ^ k) = - 4tT X k 2 , T (mag)£ (k, k) = 0. (9) 

As one see, the electric dipole gives equal contributions into transversal and longitudinal components 
of the scattering operator in momentum representation for forward scattering direction but the 
magnetic dipole contributes into the transversal component only. As consequence the difference, 
M t (k) - M e (k) = - Anfixk 2 , includes contribution of magnetic dipole scattering only. As result the 
Eq.([6]) for effective magnetic permeability takes a form 

1 47rf 1 Y , , 

1 = — ? 10 

Mf [1 - (4tt/3)/i77] 2 



7 Effective diamagnetism 

The electric and magnetic susceptibilities of spherical particle with radius tq small compared to 
the wavelength in the host medium are given according to [131 122] by 77 = Tq and x = ~ (l/2) r o> 
respectively, in the limit of perfectly reflection because of high dielectric permittivity or conductivity 
of the particle. These electric and magnetic susceptibilities of a spherical particle lead to effective 
diamagnetism of random medium consisting of such particles due to the Eq. ffTOj) . In the limit of 
small filling fraction, fxr^ « 1, these equation gives formula /i e // wl + Anfix, obtained in [3] yet. 
The Eq.( [10l . with ignoring contribution of electric susceptibility in the right hand side denominator, 
coincides in the form with a formula [23] for electron diamagnetism effect. The denominator in 
the right hand side of the Eq.( fT0l) leads to the diamagnetism enhancement because existence of 
effective dielectric permittivity given according to Eq.(|5]) by (ef// — £ o) I ( £ lff + ^ £ o) — (4tt)/3) 
fi r). One can say that diamagnetism is enhanced by appearance of additional effective dielectric 
displacement current in the medium. 
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8 Conclusion 



In this paper the physical nature of diamagnetic property in the effective magnetic permeability at 
coherent electromagnetic wave multiple scattering by statistical ensemble of independent perfectly 
reflected non-magnetic small spherical particles has been revealed. The diamagnetic physical na- 
ture consists in that the electric dipole scattering by a particle gives equal contributions into 
transversal and longitudinal components of the effective permittivity whereas the magnetic dipole 
scattering contributes into the transversal component only. This observation is made with the aid 
of generalized Lerentz-Lorenz formula technique for effective dielectric permittivity with spatial 
dispersion. 
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